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We study the spin-orbit coupling in metallic carbon nanotubes (CNTs) within the many-body 
Tomonaga-Luttinger liquid (TLL) framework. For a well defined sub-class of metallic CNTs, that 
contains both achiral zig-zag as well as a sub-set of chiral tubes, an effective low energy field 
theory description is derived. We aim to describe systems at finite dopings, but close to the charge 
neutrality point (commensurability). A new regime is identified where the spin-orbit coupling leads 
to an inverted hierarchy of mini-gaps of bosonic modes. We then add a proximity coupling to a 
superconducting (SC) substrate and show that the only order parameter that is supported within 
the novel, spin-orbit induced phase is a topologically trivial s-SC. 


I. INTRODUCTION 

In the past few years, we have witnessed a renewed 
interest in superconducting (SC) proximity effects in ID 
systems. The reason why this topic is in the forefront 
of condensed matter research was the discoverj^^ that 
a SC with a topologically non-trivial order parameter is 
able to support the long sought Majorana surface stated. 
Moreover, it was showiP that the nontrivial SC can be 
artificially created by a proximity coupling of a trivial su¬ 
perconductor with a ID wire that has a substantial spin- 
orbit coupling. While the first experimental signatured 
that such a device can indeed support Majoranas fuelled 
the interest of the community, at t he same time ques¬ 
tions about the role of disordei!^, low dimensionality 
breaking and electron-electron interactiond were raised. 
To avoid at least the first two issues one may consider 
a carbon nanotube (CNT), a self organized, strictly ID 
system that nowadays can be produced with ultra-clean 
quality. However there is still an issue of interactions 
and moreover one could wonder if the peculiar spin-orbit 
coupling, that is present in CNT, can produce topolog¬ 
ically non-trivial proximity effect. The answer to this 
questions turns out to be negative and this is one of the 
main results of this paper. 

The price of moving from a simple wire, with e.g. 
cubic structure, to a CNT is that, in the latter case, 
one deals with a highly non-trivial mapping between real 
and reciprocal space structures. The low energy physics 
of a nanotube can be derived from that of a hexago¬ 
nal graphene lattice by imposing a quantization condi¬ 
tion along the CNT circumference. For concreteness, we 
consider a CNT with a chiral vector (n, m) such that 
(u —to) mod 3 = 0. Then, within the sub-bands that fol¬ 
low from circumferential momentum quantization, there 
exist a sub-band which falls very close to the Dirac points 
K, K' of a graphene reciprocal space. The nanotube is 
metallic and the vicinities of the two distinct Dirac points 
are called valleys. More refined analysis includes a cur¬ 
vature induced shiflP away from Dirac points Acurv as 
well as a spin-orbit couplin^i^ that, in the sublattice 
basis, have both diagonal A'qq and non-diagonal Aso 
component^ii^. The spin-orbit coupling is a subject of 


particular interest due to its peculiar nature, with larger 
non-diagonal Aso component. 

It is tempting to incorporate the spin-orbit couplings 
(and Acurv) on a single-particle level because then their 
only effe ct is to change the band structure. So far, all 
attempt j^^F3 | .j-g address non-trivial proximity effects in 
CNT were based on such single particle framework. How¬ 
ever, neglecting the electron-electron interactions V (q) 
would have been justified only if they were a tiny pertur¬ 
bation added on the top of Aso and Acurv In re ality: 
U(g « 0) ~ 0.3eU and Vlg^2\K\) ~ while 

Acurv ~ Aso < lTOeV^ ^^b6 | i7 | gg gj^g faces ex¬ 
actly the opposite hierarchy of energy scales. Also, at a 
more fundamental level, a key property of ID systems is 
that even upon introducing an infinitesimally small V(q), 
their low energy description must be given in terms of 
collective excitationJi^. A carbon nanotube (CNT) is no 
exception from this general principle. A well established 
fact is that t he ve locity of charge fluctuations is strongly 
renormalizecP^EH. This is one manifestation of strong 
correlations in the physics of CNTs and it implies that 
a naive refermionization back to the original electrons’ 
framework is not allowed. 

It is then an important task to incorporate the effects 
of spin-orbit coupling into a proper many body descrip¬ 
tion of CNT. To this end a few partial problems have al¬ 
ready been solved. In Ref llOl under an assumption that 
there exists a mini-gap in the single particle spectrum, it 
has been shown that the diagonal component ^SO is able 
to shift velocities and TLL parameters of all TLL modes. 
This shift can be understood (see discussion of Eq|^ if 
one remembers that the on-site component is uniform in 
space, thus it has a density-density form. Furthermore, a 
detailed analysis of Acurv term (and interaction induced 
terms of the same form) done for a zig-zag tube, exactly 
at half filling, was done in RefljQl A crucial assump¬ 
tion was that the system is deep inside a Mott insulating 
phase. The aim our work is to go beyond this special 
case and study a new physics generated by the A^o away 
from commensurability. 

A further novelty is that a sub-set of chiral tubes is 
also covered. Apart from extending the range of validity, 
this also erases any constraints between Acurv and A^o- 
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For instance Acurv can be varied by a tube’s twistP^ (not 
possible for achiral CNTs) or, due to absence of a lattice 
inversion symmetry, an unprotected Ago can be modi¬ 
fied by higher order scattering processes. This versatility 
allows us to freely tune the parameters of our model. 

The paper is organized as follows. In Sec|TT]we identify 
a class of chiral tubes where our theory applies and then 
term by term we introduce a description within the ID 
framework. The following section Sec |III| aims to derive 
an effective low energy description in a renormalization 
group (RG) spirit. It is divided in two parts, a high 
energy part Sec |III A which is dominated by holon be¬ 
haviour, while in Sec lllf^ dedicated to the lower energies 
we use adiabatic approximation and focus on gap open¬ 
ing in the spin/valley modes. Then, in Sec |IV[ we check 
the influence of spectral gaps on superconducting prox¬ 
imity effects. Finally, in SecjVj we discuss an issue of 
experimental detection of the gaps, an influence of the 
other symmetry breaking terms, like e.g. valley-mixing 
term, and other SC orderings proposed for nanotubes. 
The paper is closed with conclusions, SecjVTj and two ap¬ 
pendices that contains estimates for a holon expectation 
value and for a proximity hybridization with a substrate. 


II. CNT AS A TWO LEG LADDER 

The hamiltonian of a CNT can be written as: 

H = Ho + H, + Hd^ + Hso ( 1 ) 

where Hq is a TLL hamiltonian (see Eq§, H, (Eq§) 
contains the many-body interactions with large momen¬ 
tum exchange, (Eqj^ is a dimerization potential in¬ 
troduced to capture Acurv, in a way following Ref Mil 
and the last, new term H^o (Eqj^ contains the Aso- 
The real space hamiltonian in fermionic language 
reads: 

Ho+di+so = X! (2) 

^ d 

where we have taken a nearest neighbour hopping on a 
bipartite lattice. The at operator creates an electron on 
the lattice site A with coordinates f. The summations go 
over all lattice sites positions r and all nearest neighbours 
d, thus d is a linear combination of hexagonal lattice vec¬ 
tors. Due to curvature effects the hopping parameter, a 
complex number tZ is anisotropic in r-space and spin 
dependent. On the top of Eqj^ one adds an electron- 
electron interaction which has a Coulomb character. In 
order to extract the low energy physics one turns to the 
reciprocal space description with a momenta , ky di¬ 
rected along CNT’s axis and circumference respectively. 
For a chiral tube both these axis make a finite angle with 
a helical line along which the graphene lattice is folded. 
The resulting band structure is illustrated in Figj^ The 
two cones, that are characteristic of hexagonal lattice, 



FIG. 1. Low energy band structnre of a CNT that falls within 
the zig-zag (like) class. As pointed out in Ref 1 141 due to 
anisotropy of tZ ^ the cones are slightly shifted away from 
K, K' points of a reciprocal nnit cell. The plane cross sections 
are due to circumferential quantization condition, only those 
valnes of quantized ky that are closest to Dirac points are 
shown. On each cone there are two dispersions E{kx). The 
mean shift away from Dirac point opens a gap Acurv, while 
a split between the two dispersions is proportional to Aso. 
Small tilt of the planes is proportional to A!,o, in R.ef llll details 
concerning these spin-orbit effects are given. 


are cut in slices that stem from the circumferential quan¬ 
tization set on ky. 

We put a chemical potential close enough to the Dirac 
points such that in the following we can restrict ourself 
only to the lowest lying sub-bands in each of the two val¬ 
leys. Creation operator are assigned to these states, 
where ct is a spin index, an index a = K,K' and fc is a 
component of an electron momentum along a CNT, thus 
ID physics is implicitly assumed. Then two Fermi points 
are present near each Dirac point and this leads to a sys¬ 
tem with overall four Fermi points. It must be the two 
leg ladder model that describes the low energy physics 
for this band structure. An exact mapping between real 
space and ct has been found for achiral armchaiil^ and 
zig-zag tubeP^. We take a closer look at the later ones as 
these can accommodate finite Acurv and Ago, the subject 
of this study. The zig-zag CNT is mapped onto a ladder 
with an interchain t_L = 0 and this allows us to identify 
chains (of an abstract ladder) with valleys (of graphene). 

The validity of this simple mapping can be extended 
also onto a sub-class of chiral tubes. In a recent workP^ 
we have showed that it is possible to distinguish a class of 
tubes defined by a condition (n — m) / gdc{m, n) mod 3 A 
0, that have two pairs of Fermi points located around 
K_i A 0) K\\ ~ ttiEit is similar to the zig-zag 
CNT. In RefOni we considered an infinitely sharp, lo¬ 
cal chemical potential, an extra term in the hamiltonian 
~ — xop(x) with po ^ oo and p(x) is an electronic 
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density, a Fourier transform of For the zig¬ 

zag like tubes a response to such potential is a reflection 
matrix that is strictly diagonal in the valley space. From 
this it follows that an operation {x = a;o)|'kfe) (where 
|d>fc) is an eigenstate) is diagonal in the valley space. We 
apply creation operation infinitely many times along a 
CNT to find that J dx is also diagonal in val¬ 

ley space which implies that a valley=chain description, 
with t± = 0 , should be valid for these chiral tubes, at 
least in the regime close to the Dirac points (fc|| « 0) 
which is of interest in this study. To quantify the cri¬ 
terion, by analogy with commensurate-incommensurate 
transitiorP^, we notice a competition between the inter¬ 
chain interaction terms ~ gic cos (j)p- (see below for defi¬ 
nition of bosonic fields and Eqj^ and the inter-chain hop¬ 
ping (present for fc|| 7 ^ 0 ) that upon bozonization gives a 
term ^ t± cos 0p_. These bosonic expressions, that con¬ 
tain two canonically conjugated fields, suggest that the 
following criterion for /c|| can be given tsin(fc||a) < gi. 
Substituting numerical values, this implies that our rea¬ 
soning can be safely applied when the doping S < 0.03. 
The fact that a CNT can be described as valley=chain 
ladder is enough to apply to the results of this work. 

We go directly to the bosonization description of the 
lowest sub-band fermions. We follow a standard proce¬ 
dure. First, one extracts the long wavelength behaviour 
around the Fermi points: 


i’aaix) = exp{ikFx)iljR^a{x) + exp{-ikFx)il)Laa{x) 

where we have written the formula in terms of a real 
space field 'ipaaix) which is an eigenvalue of the of the 
second quantization operators Cg-a(x) (Fourier transform 
of cl g.^), in the Fock space of the coherent states. Then 
one focuses on the slow components of the fluctuations 
around Fermi points and introduces the collective bosonic 
fields: 


'^Fi.Laaix') ^R.Lcra ^ ^^P(^ [^cra (^) ^ ^(7a(^)]) (H) 

Zira 

where HR^Laa, the Klein factors, ensure proper anti¬ 


commutation relations. The collective fields can be also 
expressed directly through the real-space density oper¬ 
ator defined before (for the valley-diagonalization argu¬ 
ment), for instance 4>aa{x) = —irVpgaix). Finally one 
turns to a total/transverse basis by a transformation 

4>p± = — — 4>cr± = ^ 

Four collective modes (j)^, (and canonically conjugate 
^jy) are present: total/transverse charge/spin modes {f = 
p±, cr±). The total charge mode p+ is sometimes called a 
holon as it contains an electric charge of a hole, while the 
other three modes are neutral and contain the spin/valley 
component. With these bosonic modes defined, we can 
now write down each part of Eqjl] The Hq reads: 


Ho[^. 


E f dx 

I ^ 




( 4 ) 


The main advantage of working in the bosonization 
framework is that an entire E(g « 0 ) part of interactions 
is already included in Eq[^ Since in CNTs the inter¬ 
actions have a long range Coulomb character, the small 
momentum exchange interactions are much larger than 
those with large momentum exchange. The Coulomb in¬ 
teractions bosonize as: 


Hcou = — J dx J dx'V(r - r')dx'pp+{x)dF4>p+ix') 

Clearly, only the total charge mode (holon) is affected. 
Because of this, the holons’ velocity can be up to 
five times larger than Vf, while iFp+ can be as small as 
0.25. Velocities of all the other so called neutral modes, 
remain « Vf- 

The large momentum exchange part of electron- 
electron interactions, where V{q « \K\)., or in other 
words the part that cannot be written in a density- 
density form that is 7 ^ f dx f dx'p(x)p(x'), adds several 
non-linear terms: 




= ^ / drgs cos(2</>p+ - 2Sx)(cos(2c/>p^) + cos(2(/>cr+) + cos{2(j)a-) + cos(26»o--)) - gic cos{2(j)p_) cos{24>a+) 

2[Tray J 

- g 2 c cos(2(/)p_) cos(2(;io--) + gia cos(20cr-i-) cos(20£r-) + g\\c cos(20p_) cos(26'o--) + gi cos{2(l>^+) cos(20<,r-)) (5) 


where the backscattering terms, with spin and/or val¬ 
ley index change in the process, are indicated as gi, 2 i- We 
use notation from Refill] and convention for the Klein 
factors as in ReflTSl and The only difference is that 
the Refl24l is dedicated to two-leg ladders with large t_L 


(more customary case) while here t_L = 0 but a finite 
inter-chain interaction V± is present. To transfer be¬ 
tween the two models it is enough to make an interchange 
cos 2(j)p- cos 20p- in Eqj^ Terms ~ g^ in Eqj^ are 

umklapp scattering terms which transfer two left movers 
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into right movers (or vice-versa). This requires commen- 
surability with the lattice, obeyed at half filling, while at 
finite doping S these are gradually suppressed. 

Additional terms, dimerization and spin-orbit cou¬ 
pling, are present because the C 3 symmetry of the un¬ 
derlying graphene lattice is broken upon wrapping. A 
cr* — TT* hybridization, induced by wrapping, changes the 
hopping amplitude along the tube circumference and this 
shifts the position of the Dirac pointfl^. The lowest en¬ 
ergy sub-bands are defined independently by the quan¬ 
tization condition along the tube circumference, so they 
are now shifted with respect to the new Dirac cones. This 
effectively results in an opening of a so-called mini-gap 
in the spectrum, the Acurv As it was proven in Ref l20l 
this effect can be grasped by introducing a dimerization 
potential into the effective ID hamiltonian of CNT, Eqli] 
Such a term, the Peierls term, is well known in ID sys¬ 
tems, it is exactly solvable via Bogoliubov transformation 
in the particle-hole channel and leads to a gap opening 
in the s ing le particle spectrum. In bosonization language 
it reads^Sl 


Hm = 


/ dxgdi— cosi(j>n+—5x) cos cos (j)cy+ cos <)>„- 
— sin(0p+ — Sx) sin sin ^ 0 - 1 - sin tpa -) (6) 


the tight-binding model. As a result, in the single par¬ 
ticle picture. Ago adds a spin/valley dependent compo¬ 
nent to the mini-gaps, see Figj^ By reasoning along the 
same lines like for the curvature term, this can be inter¬ 
preted as an extra spin-valley dependent single-particle 
backscattering. The Ago term is then expected to have 
a form similar to Eq|^ with the only difference that the 
left/right mixing term now involves the z-Pauli matrices 
in spin and valley space^^: 

Oso ~ c\j^^Crk^ — 

The spin-orbit coupling is expressed in the spin-valley 
basis because of the intricate topological origin of the 
effectpSIIIIISIlTl electrons of opposite valleys are precess- 
ing along the helical lines of opposite twist. However we 
have established that, within our effective two leg ladder 
description, the valley degree of freedom can be associ¬ 
ated with the chain index. Then, in Eq|^ a = K, K' 
and thanks to that Oso has a simple bosonized expres¬ 
sion. Finally, the spin-orbit term that is off-diagonal in 
the sub-lattice space, asks to choose a bond (not an on¬ 
site) operator to be hermitian. These few constraints are 
enough to deduce the following form of spin-orbit term 
Ago in the bosonic language: 


where gd = VdilVp generates mini-gaps Acurv at K,K' 
points. From the Bogoliubov transformation, done for 
an alternating potential in a single-particle limit, we 
know that in the lowest order the relation is simple 
Vdi = Acurv However Vdi incorporates also further 
terms, the staggered potential terms that are produced 
in the course of the RG flovP^. The sole term Vdi, 
EqlU written in bosonization language, contains a sort 
of ’’frustration”: there is a competition between terms 
perfectly compensating each other. Sines and cosines 
wish to lock (pi, fields at different minima. When Vdi 
dominates the physics, then the bosonic framework is in¬ 
appropriate, instead one should turn back to the original 
fermions (to obtain the Peierls transition). But this sim¬ 
ple prescription does not work if there are other terms, 
like electron-electron interactions, present as well. Then 
it is necessary to write Hdi (and Hgo) in the bosoniza¬ 
tion language in order to take advantage of the adiabatic 
approximatiorPSEmU and separate out the influence of 
the fast (pp+ field. In RefMH the ’’frustration” problem 
was solved by considering a regime dominated by the 
umklapp scattering (deep inside the Mott phase) which 
favors cos (pp+ and then also other cosines automatically 
follow. Below we show a different mechanism that is able 
to lift the ” frustration”. 

The spin-orbit coupling shifts band dispersions away 
from the Dirac points by an amount that depends on the 
spin/valley degree of freedom of a fermion, in an opposite 
direction for electrons with opposite helicitieJ^^. Alterna¬ 
tively, this phenomenon can be seen as a spin-dependent 
variation of a mini-gap in the spectrum around the point 
where bonding and anti-bonding bands used to cross in 


i?so = J dxgso{-COs{(pp+-Sx)cOS(pp- cos pcr-\-cos pa- 

-\- sin((/p_|_ — Sx) sinpp- sinpa-h sin pa-) (7) 


One immediately notices that thanks to an opposite 
sign of the two terms in Eqj^ the gso is able to lift the 
’’frustration” present in the sole Hdi- 


III. RG TREATMENT OF COSINE TERMS 


As usual in the RG procedure, we inspect how the pa¬ 
rameters of the hamiltonian are effectively changing upon 
integrating out high-energy degrees of freedom. The RG 
flow is divided in two stages: the first when the doping is 
negligible and the system flows like if it was at commen¬ 
surate filling, the second when doping is significant and 
only the backscattering terms in Eqj^ should be kept. 


A. High energy RG flow 

The first stage of RG flow stops at energy scale A' that 
is defined by the condition (5[A'] = I. Above this energy 
RG is dominated by the umklapp and dimerization/spin- 
orbit terms whose perturbative, single loop, RG equa¬ 
tions read: 


<73 = 353(1 - Kp +) ( 8 ) 

9 d,so = gd,so{‘^ — {Kp+ + Kp- -\- Ka+ -|- Ka-)/4:) (9) 
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where, in the first equation, we used the fact that 
Kp- « Kcr+ ~ Ku- ~ 1 , otherwise three different equa¬ 
tions for three different umklapp channels would need 
to be given. The reason why Eqj^ dominate is be¬ 
cause in CNT, in the UV limit, a relevant parameter 
range is 0.2 < Kp+ 1 thus one can safely assume 
— 1 | ^ \K^^p+ — 1 | and then all terms that contain 
the (/)p+ mode are much more relevant than others. The 
umklapp has a scaling dimension da = 1 — dfp+ while the 
gd and gso are even more relevant with dd = 1.25—Kp^ /4. 


The RG flow of other non-linear terms is determined 


by the following equations: 


9ic — 5ic(2 (yKp— + /^CT+)) 

(10) 

92c = 52c(2 — {Kp- + Kfj_)) 

(11) 

gia = gia{2-{K-^+K„+)) 

(12) 

9\\c = 5||c(2 - + Kp_)) 

(13) 

g3=gi{2-{Ko-+K^+)) 

(14) 


(15) 


While this flow is much slower in the first stage of RG, in 
the second stage of RG EqfT0| becomes the driving force. 
The TLL parameters are also renormalized: 


^P+ =-2^p+^'^93 + 9d + 9so)JoiS) (16) 

^P- = ~2^p-(’^o{S){g3 + 9d + 9^o) + 9ic + 92c + +5j|c) (11") 

^cr+ = + gd + gso) + 5ic + 5ia + fli) (18) 


k,_=-^Kl_{JoiS){gj + 

where Jo{6) is a Bessel function of the first kind (we take 

The bare (initial) amplitudes of the exchange terms 
in Eq{^ are small but hnite and were thoroughly cal¬ 
culated in ReffTSl In that language: gic = <?i = /, 
52 c = b — f and gia = g\\c = b, where b, f are ampli¬ 
tudes of large momentum scattering processes computed 
on a microscopic CNT lattice for armchair tube. The 
estimate &,/ « (0.05,0.1)E(g = 0 ) « (0.005,0.01)Vf 

was given and in our chiral case we are likely to be 
close to the upper limit since in a less symmetric lat¬ 
tice certain cancellation between real space Coulomb in¬ 
teractions are not exact. On the top of it, in our non¬ 
armchair case, there is a contribution from a coupling 
between orbital momenta of two electrons. It enhances 
gid gi: 9\\c (s- ferro-orbital configuration of initial orbital 
momenta /io implies that the two carriers will repel each 
other) and reduces g 2 c,gia (an antiferro-orbital conhg- 
uration of initial orbital momenta g,o) . In CNTs go 
can be an order of magnitude largep23 than gs which 
makes this unusual contribution to electron-electron in¬ 
teractions worth considering. To estimate it we can com¬ 
pare it with Ago « lAtolld-sl ^ ImeE. A^o originates 
from similar mechanism, an interaction between go and 
gB sis a, carrier moves along a helical line of a CNT. 
The umklapp terms correspond to terms with even larger 
momentum exchange, thus their initial (UV) amplitudes 
are smaller for the Coulomb-like interactions. Moreover 
their amplitude is further suppressed by a finite doping 
and this suppression is two times faster than for the Vdt 


9so) + 9i + 92c) + 9ia + 9\\c ( 19 ) 

( 20 ) 

cut-off equal to one). 


amplitude. 

Our study is dedicated to the case of a finite doping. 
Since in the later part of the paper the SC proximity 
effects are considered, we must take a model with a non¬ 
zero conductance on an interface with a substrate, thus a 
model with a constant chemical potential. Then the dop¬ 
ing is not a constant but a renormalizable quantity that 
competes with interactions. This effect we incorporate in 
the following RG equation: 

d = 5 — {Zg^-\-g^d+9lo)Ji{b) (21) 

where Ji (6) is a Bessel function of the first kind. The 
RG flow of 6 (in the first stage of RG) can produce two 
outcomes: i) (5[/] rapidly grows and when ^ 1 then 
this first stage of RG must be stopped and 53 terms in 
Eqj^and gd (and gso) terms Eqj 6 ][^ effectively drop out 
of the problem because the integrands in EqllZl contain 
rapidly oscillating terms; ii) i5[Z] rapidly drops to zero 
then the system flows to a Mott or Peierls physics, where 
a competition between g^ and gdi (and gso) determines 
the low energy properties. Case ii) can be realized only 
when (5[A] < gf which for CNTs translates into extremely 
small doping levels. Nevertheless, for a finite 6 [l] during 
RG, this competition persists and since the dimerization 
is less affected by doping then this phase should expand. 
Crucially, as we show later in Sec |IIIBT| the nature of 
the ’’dimerized” phase changes when the energy scale A' 
is of the same order or smaller than Ago. 
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FIG. 2. The RG flow of couplings that are violently relevant 
in the first stage of RG: the umklap (dashed lines) and dimer¬ 
ization (solid lines) terms. Different doping levels are shown 
by different color code (from top to bottom): from 5 = 0.002 
to (5 = 0.012 with intervals 0.002. One notices that in a point 
when flow stops, thus for I = A' the dimerization term is al¬ 
ways stronger at finite dopings even though we started with 
gd = Qso = 0.0001 and = 0.0007 (and Kp+ = 0.25). 


Close to commensurate filling, for a parameters range 
that is relevant for a CNT, we identify quite a broad 
regime where gd (and gso) dominate over g^ terms. We 
analyze several RG flows for initial parameters: gd = 
gso = 0.0001, 53 e (0.0001,0.001), Kp+ e (0.25,0.35) 
(these values are relevant for CNTs) and (5[Z = 0] S 
(0.001,0.012). Some examples of RG flows for differ¬ 
ent 6 , are given in Fig[^ We observe that both terms 
grows and, in a chosen range of parameters, the dimer¬ 
ization term is always the dominant one, even if one 
starts with (an overestimated) ratio 53 [A]/ 5 £;[A] ~ 10 . 
The flow stops for li G (5.5,7) which taking initial 
UV cut-off A = 1.5eR translates into an energy scale 
A' ^ 10“^ey that is comparable with the bare Ago- The 
values reached by 53 and gd (and gso) at A' are substan¬ 
tial gs < gd ^ 10 “^ (see Fig[^ but still below 10 °, 
thus gaps are not open yet. While these terms drops 
out of RG but in the following should be considered as a 
substantial perturbation. 


B. Physics at energies below A' ~ ImeF 
1 . Anti-adiabatic approximation 


We restrict ourself to H = Hq Hdi + Hso- At 
A' we re-an^ze the theory using the adiabatic 
approximatiorF^. To be precise we use an anti- 
adiabatic version of it to focus on the physics of three 
neutral modes. Following Reff^ we separate out 
the fast (j)p+ field using an auxiliary variable r]{x) = 
arctan[(sin 0 p_ sin (f>a-+ sin ipa -)/(cos (jip- cos 4 >o+ cos (pa-)] 
After shifting the field -G ((/)p+(a;) -I- 5x) -\- r]{x) 

the action is separable. Than for the fast field we obtain 


a sine-Gordon model: 

= Ho[(j)p+] - J drM[(j)i^p+] cos(^p+(x)) ( 22 ) 

where the mass term 


M[^,^p+] = {Vso{i) + Vdm\A' 


jl + ^ cos 2 (j)p cos 2 (j)^ 


can be obtained using identities: arctan(sin(a//?)) = 
ajip? A- /3^)“^ and (sin(/)p_ sinc^o--!-sini^^-)^ + 
(cos()>p_ cos(/)ct+cos(/)ct-)^ = 1 -b cos 2 (/)j cos 2 (/)i. 
While writing Eqj^we neglect terms ~ r]{x) (and higher 
powers) and derivatives ~ dtrj{x), which is justified in 
the adiabatic limit (slow rj{x)) and in the presence of 
substantial Wo (then g{x) —)■ 0 is justified). The Vso, as 
written in EqUl favors cosines’ over sines’ minima and 
thus provided Vso ~ A' we tend to a well defined limit 
5 —>■ 0 , variations of rj field are gradually suppressed. 


2. Effective hamiltonian for the slow fields 


For the slower fields we pro ceed by integrating out the 

is a sine-Gordon model 


22 


(/)p+. At energies A' the Eq, 
thus (cos(0p+(a;)))|t,j=A/ yf 0 (see Appendix for details). 
Then, upon expanding Mlcfi^p^] we arrive at an emer¬ 
gent non-linear term: 


H^^[<l)i^p+] = -g’^{cos{2(t>p-) cos{2(t>a+) 

+ cos{2(po-)[cos{2(j)o-^) + cos(2(/)p_)]) (23) 

where gd 5 d(cos((/)p+(a;))) and in the lowest approxi¬ 
mation the expectation value is proportional to the sym¬ 
metry breaking term {cos{4>p+{x))) 14 o- In the sign 

convention we use both hamiltonians Eqj^ and Eq |23| are 
minimized by the same combination of locked neutral 
fields, thus validating our mapping. The Eq[^ should 
be combined with the backscattering part of iJa;(Eq[^. 
The following perturbation to FIq emerges: 

Hs: = -gioCOs{2(j)p_) cos{2(j)o+) + 

C0s{2<j)o-)[-g2cC0s{2(j)p_) + gi cos(2(^,^+)]-b 

cos(2d,^_)[5iaCOs(2())„+) -b5||cCOs(2(^p_)] (24) 

The initial parameters for 2nd stage of RG flow (we 
take a new UV cut-off A') are determined by the val¬ 
ues obtained in the end of the 1st stage. The RG flow 
of Hx is a BKT flow with the parameters that falls 
close to the negative separatrixl^ (the SU(2) invariant 
line on the g — K plane, with the RG flowing straight 
away from the critical point). Then, in the first RG 
stage, g\l]/VF = g[A\/VF/{t — 5 [A]/) and one finds that 
5 ic[A'], 52 c[A'], 5 i[A'] ~ 10“^. However, some of the 
5 '[A'] terms are larger because they contain also 5 d[A'] 
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contribution (which is significant even when multiplied 
by (cos^p+)|<^=A/). To be precise: gicW] = gic[A'] + 
hcW] = 52 c[A'] + 5 d[A'], gi[A'] = gi[A'] - gd[A'] 
while all other terms are not affected. As for the TLL 
parameters, in the first part of the RG flow the umklapp, 
dimerization and spin-orbit terms all involve ~ cos{(f>p-). 
As a result the RG flow changes the TLL parameter 
Kp- downwards (already initially, at ^ = Aq, this term 
is shifted slightly below Kp_ = 1 by interaction^i^ as 
well as same holds for Kry+. Thus we conclude 

that for gic cos( 0 cr+) cos((/)p_) term we make a shift up¬ 
wards along the negative separatrix, that is both gic 
and 1 — Kp_ change upwards. Since close to the sep¬ 
aratrix the gap A = Aexp(—Vp/gic), the dependence is 
exponential. Taking quite a conservative estimate that 
both gic[A'] and 5 (i[A'] are of the same order we find 
that the exponent is reduced by a factor of two in com¬ 
parison with Ref llSI This leads to much enhanced gaps 
Ap__o-+ O.lmeR. Numerical (RG) calculations confirm 
this finding: g’lJJ,'] = 1 already for Z' « 2. A gap opens 
in the spectrum of the two bosonic modes (j)a+, 4>p-- The 
fields are locked at an energy minima (j)a-+ = 0 , 4 >p- = 0 . 
The gap value is equal to the mass of a soliton of the 
sine-Gordon model Mp_ ,j+ = 2^J2glcUp_/^TKp_, from 
this Mp-^a+ ~ O.lmey. The two estimates coincide. 

The RG flow of a— mode is more difficult to follow. 
In Hx (Eq[^ we find competing cos{(j)cr-) and cos(5ct-) 
terms which exactly compensate each other, also in the 
lowest energy sector when some modes acquire gaps. 
Moreover, this implies that dK^_/dl « 0, while to be¬ 
gin with = 1 and even accounting for the diago¬ 

nal spin-orbit coupliny^^, the A^q, does not move 
from the marginal value = 1. Thus we conclude 

that this mode is in a self-dual point, at least within 
the manifold of interaction terms we decided to take into 
account. Usually, such a situation is treated by employ¬ 
ing re-fermionizatiorP^, then separating real/imaginary 
parts as Majorana fermions, e.g. ^o,i = Re[ex:p(j)rj-(xi)], 
and finally using fusion rules to map the hamiltonian onto 
a doublet of quantum Ising chains HQ[a—] + Hx[<J—] = 
i^i i + ^Th order/disorder operators 
o’o,i;Mo,i defined as In Ref.’jTSl and E5I 

the procedure was used in the context of GNTs. The 
self dual point is equivalent, in the Ising model language, 
to CTi chain passing through criticality. The other Ising 
chain is always gapped and, by accounting for a negative 
sign of the mass term, we deduce that the order Ising 
operators (Tq have a finite amplitude, which means that 
sin (j) = 0 and sin 0 = 0 while both respective cosines are 
non-zero. This does not allow us to identify the unique 
ground state, but only to narrow down the possibilities. 
Since Kp+ < 1 it shall be DW ordering, either intra¬ 
valley GDW or intra-valley SDWz, with either bond or 
on-site character. One must remember that there are 
other ordering possibilities e.g. squared order parame¬ 
ters, with higher periodicities, which may be dominant 
when Kp^ < 0.25. Furthermore, since self-duality is not 
protected by any symmetry, one cannot exclude that due 


to some extremely tiny perturbation, not accounted in 
our generic model, a gap in actually opens. How¬ 
ever this depends on the finer details of a GNT under 
consideration and describe physics that takes place at 
energies ~ 10“®eU or belowl^, so we refrain from its fur¬ 
ther analysis here. 

Even though the exact ground state remains elusive, 
the larger gaps Mp_,M^+ that certainly open, provide 
sufficient conditions to determine the allowed proximity 
effect. 


IV. PROXIMITY EFFECTS 

The inverted hierarchy of gaps plays an important role 
in the proximity effect. This is because usually the cou¬ 
pling with the substrate and the superconducting gap 
(on the surface) are smaller than Mp-^n+. In the ap¬ 
pendix we give a brief description of the hybridization, 
in the fermionic language. To understand how these mi¬ 
croscopic considerations are linked with many body TLL 
theory, one must sum over all sites of the GNT within a 
unit cell, turn to collective fields and then express the re¬ 
sult in the two-leg ladder basis. This is a well established 
procedure, we follow RefESl to find that the singlet SG 
order operators in a zig-zag (like) GNT are: 

of*" ^ exp(z0p+)(cos (j)p- cos (l)a+ cos 0^- + *(sin <->■ cos)) 

(25) 

~ exp(i0p+)(cos())p_ sin((io.+ sin0o— -I- i(sin -o- cos)) 

(26) 

The first one, corresponds to a purely local tun¬ 
neling process (a pair is created on one site) and thus it 
is more likely to occur in type-II superconductor (short 
coherence length), than when a pair is created non- 
locally (with different phase on adjacent sites). Both 
have a scaling dimension cZa = 2 — (3 -I- R'“^()/4, 
thus they are relevant for Kp+ > 0.2. This holds when we 
assume Ki, = 1 ior u ^ p+, accounting for the fact that 
actually (in the low energy limit) Kp- < I changes the 
condition to Kp+ >0.25. It is likely that the condition 
Kp+ > 0.25 is fulfilled when a GNT lies on a conducting 
substrate which provides the screening for Goulomb in¬ 
teractions and thus reduce their range. The relevance of 
Og^ does not matter if one is deep inside the Mott phase 
and a large gap in the 4 >p+ field causes strong fluctuations 
of 0p+, thus suppressing any SC proximity effect. This 
would be the case in a system described in Ref 1201 where 
the dimerization term Hdi was governed by the Mott gap. 
In this work we have found another mechanism where the 
field ri{x) is locked by the symmetry breaking Hgo^ ^nd 
thanks to that the field 0p+ is not randomly fluctuating 
at the lowest energies. Considering the relevance of the 
SC-proximity now makes sense. 

The presence of Mp_, 117^+, or to be more precise the 
field configuration they impose, sets a constraint on the 
allowed proximity effect. If we disregard the a— mode 



for a moment, then we find that there exists one SC or¬ 
der parameter which is compatible with the locked helds 
{4>p-) = 0 and {4ia+) = 0. It is the s — SC that is 
also the most likely candidate from the microscopic view¬ 
point. This order has a topologically trivial character. 
The other Om order parameter is suppressed because 
it requires to lock the (j)^^ field at the other minimum: 
(j)u^ = tt/2. The triplet order parameters are exponen¬ 
tially suppressed because they involve the 0 ^- 1 - held which 
is canonically conjugate to the locked (l>a-+- 
As for the a— mode there are two options: 

• the mode stays on a self-duality point: Then 
cos 9cr- has a finite expectation value and s-SC is 
allowed; 

• ultimately the gap opens, at much reduced 
energies: this will be most likely a gap in the (pa- 
field, rua-- It could suppress the s-SC proximity 
effect at the lowest energies. One way to overcome 
it is to take a sufficiently large amplitude of the 
proximity induced gapAf*^ > ma-- Thanks to a 
huge difference of energy scales between the differ¬ 
ent masses the ID character of the system will be 
still protected by Mp_. The induced transition to 
the s-SC state shall have the Ising charactei!^ (one 
of Ising disordered operators p,i acquires a finite 
value at the cost of the Ising ordered operator cti ) 

In either case only the topologically trivial s-SC is al¬ 
lowed. The Majorana surface states are never allowed to 
occur. 


V. DISCUSSION 

A. Size of the spectral gap and the means of its 
detection 


The estimate for spectral gaps that we have given in 
Sec IIIB| is rather conservative, valid for a CNT embed¬ 
ded in a good dielectric, for instance a CNT suspended 
in vacuum. By introducing an extra screening, for exam¬ 
ple by placing a tube on a superconducting substrate or 
within a multi-wall CNT, one makes electron-electron in¬ 
teractions more local. In reciprocal space this increases 
the large momentum exchange component of electron- 
electron interactions, V{q ~ 2K). Then the bare am¬ 
plitudes of the backscattering terms gi in Eqj^ can grow 
substantially. Moreover, as we indicate in the context 
of the proximity effect, placing the tube on an appropri¬ 
ately chosen substrate may introduce additional periodic 
potentials that cause backscattering and adds up with 
and gso- The magnitude of gaps depend on particular 
experimental realizations and in some circumstances it 
can be detectable already at energies ~ ImeV. 

One possibility to detect the Mp-,Ma+, is to study 
the Knight shift and relaxation rate of NMR signal. The 
temperature dependence shall be a power law but at the 


energy scale corresponding to the gap one should ob¬ 
serve a change of an exponent, such an effect was indeed 
experimentalljI22MI observed but its origin was unclear. 
In our mechanism, for instance for the Knight shift we 
predict a change from {Ka+ + Ka-)l2 to Ka-I2. More¬ 
over, a known feature of the spin-valley dependent split 
Vso is tha t it ca n be varied by applying an external mag¬ 
netic fiel d^^ l ^^l Since both Mp_ and Ma+ gd I4o, 
and the spin/valley dependent part of the split in a single¬ 
particle dispersion can be varied by a magnetic field di¬ 
rected along a tube, then an anisotropic magnetic field 
dependence of spectral gaps can be taken as a hallmark 
of their many-body origin. 


B. Relation to SC order parameters proposed for 
CNTs 


The in the same form like Eq 25 was also proposed 
by Eggei^. The fermionic expression, in the reciprocal 
space, for superconducting order parameter that we in¬ 
voked reads: 


or = r 


-kK'i + - (t^i) 


kK^^-kK'i 

and it is equivalent to an inter-chain ordering as de¬ 
rived in a seminal papei^^. In the last paper it is called 
d-SC, but this should not lead to any misunderstanding, 
since we define order parameters for real-space hexago¬ 
nal lattice, what is Od for a square ladder is not nec¬ 
essarily d-wave for other underlying crystal lattice. A 
detail description of the symmetry properties for a bi¬ 
layer graphene interface is given in B e f i8^ where a tables 

as well as the non- 


of characters for the local Eq 
local O^, the Eq, 


26 


25 


were found. In particular it was 
explicitly shown that only the may contain topolog¬ 
ically non-trivial SC order. 

Furthermore, one notices that Of is different from the 
superconducting order parameters proposed previously 
for the armchair CNTs^. This is because the band struc¬ 
ture is different: the inter-band order parameter, that 
was previously prohibited due to the conservation of A:| |, 
now is allowed because in zig-zag (like) tubes the chains 
of ladder are associated with valleys and Dirac cones are 
located at ATn = 0. Moreover, if the circumferential mo¬ 
mentum is conserved, then by requiring ki = —k 2 within 
the BCS pair, we Hnd that indeed the inter-chain (inter¬ 
valley) Of*^ is favored (see Appendix for details). More¬ 
over, from a basic symmetry argument, we know that 
the inter-valley Andreev reflection is protected (vs for 
instance disorder) by time reversal symmetry. On the 
other hand the intra-valley pairing would be protected 
by the so called symplectic symmetry, but this one is al¬ 
ready broken from the very beginning by introducing the 

Aso- 

So far we have discussed the relation between Off) 
and and other uniform SC-orders proposed before. A 
novel aspect of proximity effect, that is inevitably present 












9 


in chiral tubes, is its non-uniformity. For a chiral tube 
that is rolled along the helical line one may consider the 
skew-turn to be a built-in rotation angle vs substrate 
lattice. Since strength of bonding is related to inter¬ 
atomic distance,the de Moire pattern of the substrate- 
tube hybridization appears and the proximity effect is 
not any longer uniform but instead it becomes periodic!^ 
(see Appendix for details). Such periodic proximity ef¬ 
fect is favourable for more exotic SC orders proposecP^ 
for the two leg ladder models and known as pair density 
waves. One defines a composite order parameteil^, that 
is a product of the Of ^ and some density wave. The den¬ 
sity wave shall be defined in the intra-valley channel to 
avoid a direct competition with superconductivity. One 
advantage is that one can construct an operator Of^^y 
which, albeit less relevant, depends only on (j)u+ in the 
spin sector. This allows to avoid a potential problem if 
a field is after all locked. The SC order most likely 
retains a topologically trivial character, in the sense that 
standard procedure of Refl551 again favours Of*". There 
are many other fascinating aspects of non-uniformity that 
should stimulate research in this direction. One is that 
the SC proximity effect shall be particularly strong for a 
chemical potential for which kdM = kp. This opens an 
exciting perspective of gate tuning of SC order in CNTs. 


C. Other symmetry breaking terms; valley mixing 


For completeness we comment on other backscattering 
operators, analogous to Hdi and iJgo, that can be intro¬ 
duced into the hamiltonian of a CNT. One frequently 
proposed perturbation is an inter-valley backscattering, 
the so called Akk'- No matter what the content in the 
spin-space we choose, this operator written in bosonic 
language contains 9p_ field, a field canonically conjugate 
to all cosine terms present in Eq’s lllzl This means 
that Akk' is quickly suppressed by all other terms as 
one is moving along the RG trajectory (towards L —> oo 
in quantum dot language). Moreover, even in the case 
when an extremely strong App/ is able to dominate the 
physics, since the pairing operators in Eq’s|25|26 


con¬ 


tain cos (j)p- then none of these (including the poten¬ 
tially topologically non-trivial, non-local Off) shall be 
favored. It seems that the Appi reduces the propensity 
of the system to any standard proximity effect. 


One can also ask an opposite question; what symmetry 
breaking term could potentially support the 0^ order? 
A brief inspection of all order parameters reveals that 
this is a rather exotic spiral electric field acting opposite 
on two valleys (but valley-diagonal), however this would 
need to be taken together with attractive V{q « 2\K\) 
interactions. 


VI. CONCLUSIONS 

We have shown that, for a chosen sub-set of CNTs, the 
presence of spin-orbit coupling A^o leads to a gap open¬ 
ing in the spectrum of two bosonic modes (j)p- and 4 )^+ ■ 
This drastically reduces the sub-set of proximity effects 
allowed at the lowest energies: we find that only a phase 
with a trivial topology is allowed. This statement is quite 
general as it should remain valid also upon increasing 
the interaction strength, doping, hybridization with the 
substrate and upon adding another symmetry breaking 
term App'. An extra motivation, and a broader per¬ 
spective, for this work comes from the recently synthe¬ 
sized 2D analogs of graphene: silicene, germacene and 
stanene. This gives a hope for a new class of nanotubes 
that shall be built out of atoms heavier than carbon. 
Since Aso ~ ^so (where Ago is an atomic spin-orbit cou¬ 
pling constanflAiJ) , then the fine effects predicted here can 
become orders of magnitude larger. 


Appendix A: Estimate for (cos(^p+)|ij=A' 

The dynamics of the fast field (j)p+{x) for energies ^ A' 
is rather complicated. Usually i5 ^ 1 implies that the gd 
and gso terms rapidly drop out of the problem (and ga 
as well, but by writing Eq|^we had already neglected 
gs). However, already in the simplest single mode sine 
Gordon model, the issue of how precisely the expectation 
values disappear when S becomes substantial, has proven 
to be quite-nontrivial and depends on how precisely RG 
procedure is set ufP^. Our model is much more compli¬ 
cated as the dynamic coupling with three other modes 
is present. For instance in the argument of the cosine 
one can clearly see the competition between ri{x) and Sx. 
Moreover, the amplitude of the cosine, that is M[(j)i^pp], 
shall have an extra increase when the two neutral fields 
order. 

To tackle the problem let us assume that, to begin 
with when r]{x) « 1, r]{x) dominates. At A' energy 
scale the neutral fields still fluctuate, with a velocity 
that is irrational with the holon velocity, thus M[(pi^pp] 
and the g{x) field can be considered as amplitude and 
phase of a complex random variable. Then Eq |22| can 
be interpreted as a model of a random backscattering 
in a TLL. The suppression of gd is delayed by the fact 
that for the disorder problem the scaling dimension is 
even larger ddis = 3 — 2Kpp. One can say that a 
strong enough disorder freezes the correlation function 
for I Ks A'. Such a phenomenon occurs also for incommen¬ 
surate (not necessarily random) potentials when modes 
of different velocities couple. It is then known as Aubry- 
Andre transitioiP^. The crossover is quite complex, but 
most likely as the energy scale I decreases during the 
RG, then g ^ 0 (the fluctuations cease below the energy 
~ Mp-), the randomness disappears and g 3 ,d,so[l] resume 
their flow to zero, driven by a finite doping. However, for 
energies around A', the sine-Gordon model, Eqj^ with a 
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finite, energy independent amplitude of the cosine term, 
gives a correct description. Then one can attempt to 
compute (cos0p_|_)|tj=A', in a limit when the shift goes 
to zero, by using results known from the Ising model 
in the renormalized classical regime. One may either use 
the zero temperature result ^ Ko{gicT'), where Kq is the 
modified Bessel function of the second kind and a charac¬ 
teristic time-scale is set as r' ~ 1/A', or a finite tempera¬ 
ture result BeffMl where (cos (/p+) |(.j=a' is proportional to 

erfc{\jTl^gic) with a characteristic temperature taken 

to be ksT = A' (and er/c is a complementary error func¬ 
tion). In both estimates we get (cos^p_|_)|^=A' ~ 10“^. 
We consider it as an upper limit for (cos ^p+)|t,j=A' and 
in all further calculations we take a more conservative 
value (cos ^p+A' = 10“^. 


Appendix B: Details of an overlap with a substrate 

In RefHni it is shown that every site which is in 
touch with a superconductor, upon integrating out 
the BCS condensate, acquires an emergent pairing 
potential:/ dxdyAsc{x,y){cl.{x,y)cLg.{x,y) + h.c.) where 
{x, y) are the site coordinates (interface has 2D charac¬ 
ter) and Asc{x,y) ^ t” {x^y^/Vp is a pairing strength, 
with t” (x) a hybridization between CNT and a substrate. 
Let us consider a process of creation of a Cooper pair in¬ 
side a CNT: . For a moment we need to take a 2D 

All Ai2 

ki because we keep interfaces’ 2D character. When the 
pair is created in two different valleys (an inter-valley 
term), it is compatible with the standard s-wave BCS 
pairs in the substrate where ki = — ^ 2 . Contrarily, the 
intra-valley term does not conserve momenta since then 
ki± = k 2 ± ± 2K±. Thus this second process will be 
suppressed when ki± is a conserved quantity during the 


tunneling process. We can try to quantify the condi¬ 
tion for conservation of the circumferential momenta. We 
take the hybridization t” [x, y) to be a Gaussian with a 
width proportional to the nanotube radius: 6y = ai,R, 
where at is some proportionality constant and ?/ is a di¬ 
rection along tube’s circumference. This relation sim¬ 
ply encodes the fact that for broader tubes there are 
more carbon atoms that can build a covalent bond with 
a substrate. Finite 5y produces a momentum resolution 
Skperp ~ l/<5y. The two valleys can be distinguished pro¬ 
vided the real-space Gaussian is broad enough, that is 
Skperp < \K\ R > a/\K\. When this condition is ful¬ 
filled one can consider valley index and thus k± to be 
a conserved quantity in a substrate-CNT tunneling pro¬ 
cess. 

The microscopic model also allows us to take a closer 
look at the non-uniformity of the t”(a:). For a chiral 
tube the hexagonal lattice makes consecutive skew-turns 
around the central axis of the tube. Then looking from 
the top it is very much like a sequence of tilted hexagons 
{Sy > with a graphene lattice constant, is as¬ 

sumed). If one puts two hexagonal lattices one on the 
top of another and rotate (or re-scale) one of them then 
he obtains the periodic de Moire pattern. Re-scaling is 
necessary only when the substrate is a crystal different 
from graphene. We conclude that a chiral CNT placed 
on the top of a 2D surface, gives an effective hybridiza¬ 
tion r (x) = dyt” (x, y) that is not constant along the 
tube but varies and these variations are the strongest for 
smaller tubes where the effect is not averaged out by in¬ 
tegration over large Sy. For scaling factor between two 
lattices equal to one (e.g. both based on graphene) one 
find^^ that the angle between the kdM and the CNT ba¬ 
sis is 7r/2 and indeed the hybridization t"{x) along ID 
profile is periodic. The patterns’ periodicity depend on 
the chiral angle, for small chiral angles very small \kdM\ 
can be reached. 
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